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Monte Carlo Simulation and Stochastic FEA are used to predict 
randomness in the free vibration response of thin unsymmetrically lam- 
inated beams. For the present study, it is assumed that randomness in 
the response is only caused by uncertainties in the ply orientations. The 
ply orientations may become random or uncertain during the manufac- 
turing process. A new 16-dof beam element, based on the first-order 
shear deformation beam theory, is used to study the stochastic nature of 
the natural frequencies. Using variational principles, the element stiff- 
ness matrix and mass matrix are obtained through analytical integration. 

Using a random sequence a large data set is generated, containing pos- 
sible random ply-orientations. This data is assumed to be symmetric. 

The stochastic- based finite element model for free vibrations predicts 
the relation between the randomness in fundamental natural frequen- 
cies and the randomness in ply-orientation. The sensitivity derivatives 
are calculated numerically through an exact formulation. The squared 
fundamental natural frequencies are expressed in terms of deterministic 
and probabilistic quantities, allowing to determine how sensitive they 
are to variations in ply angles. The predicted mean- valued fundamental 
natural frequency squared and the variance of the present model are in 
good agreement with Monte Carlo Simulation. Results, also, show that 
variations between ±5° in ply-angles can affect free vibration response of 
unsymmetrically and symmetrically laminated beams. 


INTRODUCTION 

In recent years, there has been an increasing de- 
mand for laminated composite materials in aircraft 
structures. The main reasons are that the compos- 
ites posses the following characteristics: lightweight, 
cost-effective, and can handle different strengths in 
different directions 1 . However, these materials offer 
quite a few challenges to structural engineers. Be- 
cause of their inherent complexity, laminated struc- 
tures can be difficult to manufacture according to 
their exact design specifications, resulting in un- 
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wanted uncertainties. 

The design and analysis using conventional mate- 
rials is easier than those using composites because, 
for conventional materials, both material and most 
geometric properties have either little or well known 
variation from their nominal value. On the other 
hand, the same cannot be said for the design of struc- 
tures using laminated composite materials. The 
understanding of uncertainties in the laminate com- 
posite structures is highly important for an accurate 
design and analysis of aerospace and other structures 
using composite materials. 

These uncertainties are defined as random- 
ness from non-cognitive sources involving physical 
stochastic likelihood and human factors”. This ran- 
domness can occur in each layer and involve quan- 
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tities such as: ply-orientations, thickness, density, 
and material properties among others. Such varia- 
tions can affect the behavior of the structure. 

What has intrigued many engineers to study 
structures with laminated composites is the com- 
plexity of these materials. Since it is costly to an- 
alyze a composite structure as a three-dimensional 
solid, analysis of many composite structures can be 
performed using laminated one and two dimensional 
theories, such as beam theory. A review of vari- 
ous available theories for analyzing laminated beams 
is given by Raciti and Kapania 2 3 . Earlier, a 12-dof 
element was developed and formulated for deter- 
ministic symmetric laminated beams to study their 
static and dynamic behaviors 4 . Later a 20-dof el- 
ement (Kapania-Raciti Element) was developed to 
study static, free vibration, buckling, and nonlinear 
vibrational analysis of unsymmetrically laminated 
beams 5 ’ 6 . In both works the effects of uncertain- 
ties were not considered. 

In considering uncertainties in a system, one can 
consider three different approaches: (i) probabilistic 
methods, (ii) fuzzy set or possibility-based meth- 
ods, and (iii) antioptimization 7 . Here, however, only 
the probabilistic approach is considered. The prob- 
abilistic analysis can be performed using either an 
analytical or a computational approach. An ana- 
lytical approach would be most accurate although 
cumbersome and impractical except for very sim- 
ple systems. However, with the availability of ex- 
tremely fast computers, the finite element method 
has become a widely used technique for design and 
analysis in engineering. Therefore, a finite element 
approach using Monte Carlo Simulation is developed 
to take into account the stochastic nature of the ply- 
orientations. 

Work done by Vinckenroy 8 presents a new tech- 
nique to analyze these structures by combining the 
stochastic analysis and the finite element method in 
structural design. Agrawai et al. 9 used a wavelet- 
based stochastic analysis to analyze isotropic beam 
structures. 

Librescu et al. 10 studied the free vibration and 
reliability of cantilever composite beams featuring 
structural uncertainties. They used a Stochastic 
Rayleigh-Ritz formulation. However, to the best of 
the authors’ knowledge, no work has been found re- 
garding the effect of uncertainties incurring in the 
ply-orientations on the natural frequencies of thin 
unsymmetrically laminated beams using Stochastic 
Finite Element Analysis. 

The primary goal of this investigation is to pre- 
dict the dynamic behavior of thin unsymmetrically 


laminated beams with uncertainties. In order to 
study the stochastic nature of the dynamic response 
of such beams, a new 16-dof element is developed 
using first-order shear deformation beam theory to 
account for uncertainties. Only an overview of this 
element is presented here, and a more detailed anal- 
ysis of this new 16-dof element will be presented 
elsewhere. Using this element the free vibration 
analysis is performed for only those uncertainties as- 
sociated with layer-wise ply-orientations. 

16-DOF LAMINATED BEAM ELEMENT 
An overview 

The present 16-dof laminated beam element takes 
into account the existence of various coupling ef- 
fects, which play a major roll in laminate composite 
materials. This element is valid for the analysis 
of both symmetric and unsymmetrically laminated 
beams. The motivation to develop a new beam ele- 
ment was to have a formulation consistent with the 
first-order shear deformation beam theory (FSDT), 
able to analyze unsymmetrically laminated beams, 
and that would account for most of the uncertainties 
involved in a thin laminated beam when modelled 
using FSDT. In the present work, the reference sys- 
tem of coordinate is such that the x - axis lies 
along the length of the beam and the z - axis is 
placed at the mid-surface measuring the transverse 
displacements. The present work assumes that the 
x - z plane divides the beam in two identical parts: 
in other words, material, geometry, and loading are 
symmetric about the x - 2 plane. When considering 
twisting and ignoring in-plane shear, the displace- 
ment field for the first-order shear deformation beam 
theory, in the defined reference system, can be ex- 
pressed as follows 


U{x,y,z) 

= u(x) -r zo(x) 

(la) 

V{x,y,z) 

= 0 

(lb) 

W(x,y,z) 

= tx(x) - yr{x) 

(1c) 


Therefore, the finite element formulation considers 
eight degrees of freedom at each node: axial dis- 
placement u, transverse deflection u\ rotation of the 
transverse normals <£, twist angle r, and their deriva- 
tives with respect to x. These nodal displacements 
are denoted as 

{< 7 } = {u<u\w,w\<P>q'.t,t'} T ( 2 ) 

The deflection behavior of the beam element for the 
first-order theory including transverse twist effect is 
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described as follows 

u(x) — N\ Til ■+* No u[ + N$ y >2 4* N* u' 2 
w(x) = N\ w\ + N 2 u/j + N 3 u /2 + N 4 u)o /^\ 
(j>(x) — iVj 0! + N 2 <P\ + N$(j >2 + N$(j >2 

r(x) = N\ T\ + N 2 T{ + jYj T2 + A r 4 T2 

where the shape functions N t (x) are the well known 
Hermitian polynomials. 

Constitutive Material Law 

In this investigation, all laminated composites are 
considered as orthotropic materials. Also, consider- 
ing a state of plane stress and eliminating e zz from 
the stress-strain relationship, the reduced material 
coefficients are expressed as Q ty Since in laminated 
composites each ply may have different orientation, 
the stresses are expressed in terms of an arbitrary 
angle 0 and the transformed stress-strain relation- 
ship takes the form of 11 


is identically zero and the coupling between bending 
and stretching vanish. However, for unsymmetri- 
cally laminated beams the coupling cannot be ig- 
nored and it must be included in the analysis. In the 
presence of uncertainties, laminate composite struc- 
tures are no longer symmetric and the analysis of 
unsymmetricaliy laminated structures is a more ac- 
curate one. 

Laminate Constitutive Relations 

The basic constitutive relation is 

{N} = [D\ {e} (9) 

where {N} T is the stress resultant vector, [D ] is the 
bending-stiffness matrix, and {e} 7 is the stretching 
and bending strain vector. 

The displacement field, Eq. (1), suggests that 

2e ly = 7xy-t-2«xy=° U°) 



Q n Q 12 Ql6 

Q 12 Q22 ^?26 

QlG Q 26 ^66 
0 0 0 
0 0 0 


(4) 


0 b 1 { e xx 

0 0 e yy 

0 _0 < 2 e ly 

Q44 Q45 “ e y* 

Q 45 Qss J l 2eiz 


Therefore, N ly and M ly are not considered because 
they vanish in the strain energy formulation. In ad- 
dition, the present formulation assumes: 

M yy = iV yy = 0 (11) 

This leads to the following bending-stiffness matrix: 


[D 1 = 


[Du] [ Dui \ 
[Duj ] [Du, 11 } 


( 12 ) 


where 


where Q t) are the transformed stresses 1 . 

[A], [B], [D] Matrices 

The extensional matrix [A], the extensional- 
bending coupling matrix [B] , and the bending stiff- 
ness matrix [D] are of great importance in the 
present work and these are calculated as follows 11 


J^lom 

,(z k +i-Zk) = 1 . 2,6 

k— 1 

Starr. _ 

Atj = K Y, Qm (*»+« " l 'J = 3 4 > 5 

k~ 1 



(5) 

( 6 ) 

(7) 

(8) 


where K = |, the shear correction factor 12 , and 
:V/am is the total number of plies considered. When 
considering symmetrically laminated composites, [B] 


[£>/,/] 
[Du, 1 
[D,u 1 
[Du,,, | 


An 

B n 

B n 

D n 

A\2 

B \ 2 

B 12 

D 12 

A i2 

B\2 

B 12 

D\2 

A 22 

B 22 

B 22 

Do’2 


(13) 

(14) 

(15) 

(16) 


The reduced form of the bending-stiffness matrix is 
calculated as 

[Dr 1 = [Du] - (D/.//I [D,ui\~ l [ D,u 1 (17) 

Using the above expression, an equivalent 
bending-stiffness matrix [D c ] for a thin unsymmet- 
rically laminated beam is found: 

N xx ) \ D cu D c „ 0 0 

A/ IX I __ D c 13 &C22 9 0 

Q x [ 0 0 D Cyj D Cm 

Qy I 0 0 D cu D C44 

( 18 ) 



3 OF 11 


American Institute of Aeronautics and Astronautics 



where. 

The minimization of Eq. ( 21 ) results in the element 
stiffness matrix for the 16 -dof laminated beam ele- 


ment: 

A = B 22 — A 22 B22 


D Cu =(Ao2 B\2 - 2 A12 B\2 #22-r 

[K e ] = 7 7 [ B 3d \ T [D c j [B , d i dy dx ( 22 ) 

Jo j ~ 

An B \ 2 -f A\ 2 D22 - An A22 D22 )! A 


D Ci2 ~(B\ 2 B22 + B\\ 4 * A22 B\o D\ 2 ~ 

A 12 B22 B)\2 “ A22 B\\ D22 + A 12 B\2 7 > 22 )/ A 

where b Q is the width of the beam. £ e is the length 
of the beam element, and [D c ] is the equivalent 
bending-stiffness matrix. 

D Cl 2 —{B22 D\ 1 — 2 Bi 2 B22 D\ 2 + 

A22 B\2 4 “ B\2 ^22 — A22 ^11 D22 )/ A 

Element Mass Matrix 

The kinetic energy for this 16 -dof beam element 
is 

Dc SJ ~ A44 
~ A 4 5 
£>c 44 =A ,55 

T = 1 JJJ p [t/ 2 + V 2 ± lV 2 ] dV ( 23 ) 
\/ 


Note that the availability of symbolic manipula- 
tor like MATHEMATICA®Version 4 . 0 * has made it where p is the mass density. The mass matrix is 

possible to determine the above matrix analytically. obtained by substituting Eq. ( 1 ) into Eq. ( 23 ) and 

By obtaining this matrix analytically, the CPU time then taking the first variation of the kinetic energy: 
is saved, a great help in the Monte Carlo Simulation. 

Strain-Displacement Relationship QJ QJ OT 07 . 

7 7 7 7 • i j + — dw -b + -^~ dr C 24 ) 

The strains in the above formulation are related ou ow o<p or 

to the displacements as follows: 

{ s xx ) f 1 The mass matrix coefficients are obtained as follows 




Using Eqs. ( 2 ), ( 3 ) and ( 19 ) the strain-displacement 
relation can be expressed as 

{e} = [B sd | { q } (20) 

where [£*<*] is the strain-displacement matrix. 

Element Stiffness Matrix 

Using Eqs. (9) and ( 20 ), and noting that matrix 
[ D c ] has been integrated throughout the thickness, 
the strain energy for this new 16 -dof beam element 
becomes 

U= l -{q} T JJ[B, d } T [D c }[B sd \dA {<?} ( 21 ) 


t Mathematical is a registered trademark of Wolfram Re- 
search. Inc. 


Suj : Iq / Ni(x)Nj(x )dx u,+ ( 25 ) 

Jo 


f' Ni(x)Nj(x)' 
Jo 


dx I o, 


Swj : Io f N l (x)N J (x)dx w, ( 26 ) 

Jo 

6 <t> } : [ N,(x)N } (x) dx u,+ ( 27 ) 

Jo 

h f N l {x)N ] {x)dx <Pi 

Jo 

6tj :\j x J N x {x)N } (x)dx r, ( 28 ) 
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where 


-v, 


Io = bo £ P k 

(Zk + l ~ Zk) 

(29) 

k=l 



N | am 

h = bo £ P k 

k— l 

f z k+l ~ Z k\ 

(30) 

l 2 ) 

^lom 

h = K £ P k 

k— 1 

CM 

(31) 

, 3 

r °° T J 

J ’ = 12 E p 

! (Zfc+i - Zk) 

(32) 


fc= 1 


In the present work, analytical expressions for co- 
efficients of both mass and stiffness matrices are 
obtained. Once again, this greatly saves CPU time 
for the Monte Carlo Simulation as one no longer has 
to perform numerical integration for each case. 

FREE VIBRATION ANALYSIS 


where cr 2 is the variance of the random variable: and 
/i is the mean value of the random variable. In the 
present work, the random variables are considered 
as independent and are denoted as 

x = {$ u 0 (38) 

where O^s are the ply-angles. 

Various methods exist to analyze an uncertain 
unsymmetrically laminated beam by integrating 
stochastic aspects into the finite element modelling: 
perturbation techniques, Taylor Series, and Monte 
Carlo Simulation. These techniques have been used 
in the past two decades in fields involving random- 
ness. Especially, in the last decade there has been a 
growing interest in applying these methods to better 
understand laminated composite structures. More- 
over, with the finite element technique becoming so 
popular, there has been a new interest in integrat- 
ing the stochastic nature of the structure in the finite 
element analysis 15 . 


The Hamilton’s principle is used to study the dy- 
namic nature of the structure. The principle uses 
the Lagrangian which is defined as 

£ = T -U-V (33) 


where, 

T is the total kinetic energy defined by Eq. (23) 
U is the total strain energy defined by Eq. (21) 


V is the potential of the applied loads 


The Hamilton's Principle can be represented as 13 


SH 



- 8U - 6V} dt= 0 


(34) 


The free vibration analysis is obtained by setting 
8V = 0 in Eq. (34). This leads to the free- vibration 
equations of motion for the defined finite element 
problem: 

[M] {q\ + [K\ {<?} = 0 (35) 

where [A/] is the mass matrix and [A'| is the stiffness 
matrix. Further assuming a harmonic response, the 
solution of these equations results in an eigenvalue 
problem: 

[[K]-u> 2 [M]] {<?}=() (36) 


A PROBABILISTIC APPROACH 


An overview 


The present analysis assumes that the random 
data is symmetrically distributed. Therefore, it as- 
sumes normal distribution 14 : 


/(*) 



exp 



21 


(37) 


Monte Carlo Simulation 

Monte Carlo Simulation, although computation- 
ally expensive, methods are quite versatile tech- 
niques capable of handling situations where other 
methods fail to succeed. The MCS is also used 
to verify the results obtained from other methods. 
Monte Carlo Simulation methods are based on the 
use of random numbers and probability statistics to 
investigate problems. For purposes of the present 
work a random number generator is used to generate 
possible angle- variations between -5.0° and 5.0°. 

A large sample is generated and then using PDFs 
one evaluates the probability of having such values. 
The larger the number of simulations more the con- 
fidence in the probability distribution of the results 
obtained. Therefore, for the present analysis at least 
ten thousand realizations of the uncertain beam are 
performed, increasing the accuracy of the ply-angle 
distribution fitted to the sample data. 

Stochastic Eigenvalue Analysis 

The stochastic eigenvalues problem is expressed as 

[K - AM] {<£} = 0 (39) 

where K , M, A, and <p are the stochastic stiffness 
matrix, mass matrix, eigenvalues and eigenvectors, 
respectively. 

The presence of uncertainties in ply angle orienta- 
tions results in certain randomness in the extensional 
matrix [A], bending-stretching coupling matrix [B], 
and bending matrix [D] . Since the coefficients of 
these matrices are present in the equivalent bend- 
ing stiffness matrix [D c ], Eq. (18), the matrix [D c J 
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will have certain randomness associated as well and 
these uncertainties are expressed as 

[D e \ = D c ° + D c 3 (40) 


of their mean-centered zeroth-, first-, and second- 
order rates of change with respect to the random 
variables as 
Eigenvalues 


where D c ° is the equivalent bending stiffness ma- 
trix evaluated at the ply mean values and D c 3 is 
the equivalent bending stiffness matrix taking into 
account the random nature of the problem. As a 
result, these uncertainties affect the stiffness ma- 
trix. However, the mass matrix is not affected when 
only considering those uncertainties involving ply- 
orientations. Therefore, the mass matrix M is eval- 
uated at the given deterministic values. 

In problems where uncertainties are considered, 
there exists no density function describing the ran- 
dom nature of the system. The information is lim- 
ited to only the mean values of the random vari- 
ables. In such cases, perturbations techniques are 
suggested among other existent techniques 15,16 . In 
the free vibrational analysis of the present problem, 
the random nature of the stiffness matrix, eigen- 
values, and eigenvectors are studied using a Taylor 
series expansion up to second order about the mean 
of each random variable. The present formulation 
follows the approach followed by Librescu et ai. 10 
Opposed to the work done by Librescu et al., the 
present formulation can use results provided by com- 
mercial finite element codes, greatly reducing the 
number of calculations, and calculates the sensitivity 
of the eigenvalues up to the second-order approxima- 
tion. 

Therefore, the random nature of the elastic stiff- 
ness. K , is expanded in terms of the mean-centered 
zeroth-, first-, and second-order rates of change with 
respect to the random variables as 

K(x\ , £2 Xn)= ( 41 ) 

K°+±K's t+ [ -±±K>‘ S% e } 

1=1 1=1 J = 1 

K I 

lx=(^t ,/ia, ,**«) 

dK I 
d 2 K | 

dx t dxj #*r») 

“ Mi 


K° = 

K‘ = 
= 


The eigenvalues and eigenvectors are also affected 
by uncertainties in the ply angles. The perturbed 
eigenvalues and eigenvectors are expressed in terms 


A(xi,z 2 , . . . ,x n ) = 


(42) 


1=1 


ill- 




1=1 ;=i 


A 0 

K 1 


a| 

lx=(/ii ,M2 ,- ^n ) 

ax I 

dx x ix = (^l ) 

a 2 A | 

dx x dXj lx=(/il,M2 ^n) 


Eigenvectors 


{4>(Xi,X 2 , . . . ,X n )} = 


(43) 


{^°} + } g < + 

*=i 


5 x; 


.=1 j=i 


{* 0 } 





dm 

dx t ,^n) 

a 2 W I 

dx % dxj lx=(^i,M2. -m«) 


After substituting Eqs. (42), (43), and (44) into Eq. 
(39), the stochastic eigenvalue problem is expressed 
as 


k ° + £ K ' £ ' + \ £ £ K 'l s ' £ > 

i=l i=Lj=l 

* { ^° + £ ^' £ ' + •> £ £ } 

y .=1 •=! i=\ ) 

.p + i;^.+5x;x:< ! d i m ": 

\ »=1 i=l J=1 / 

if; <«> 

[ 1=1 t=l)=l J 


The uncertainties in the random variables are in 
general small. As a consequence, in the applied per- 
turbation technique it is sufficient to only consider 
the first and second derivatives of eigenvectors and 
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eigenvalues with respect to the random variables. By 
equating the zeroth order terms of £ t in Eq. (44), 
an eigenvalue problem for the mean-valued system 
is obtained. From which, the mean-centered zeroth 
derivative eigenvalues and associated eigenvectors 
are obtained 

[K° - A°M°] {<£ 0 } = 0 (45) 

Premultiplying the first and second order terms of e x 
in Eq. (44) by {<£°} r and simplifying, the following 
expressions are obtained 

{</>°} T \k[ - A'M 0 ] {<£ 0 } = 0 (46) 

{<p°} T [if-/ - A'/M°] {tf> 0 } = 0 (47) 

This results in expressions for the mean-centered 
first and second eigenvalue derivatives: 

A , {®°} T [«:,'] {»"} 

’ ' {*°} T 

y K'] y 

{*»} T [m»] (®“: 

The advantage of this method is that the eigenvalue 
problem needs to be solved only once. The sen- 
sitivity analysis is done by using results from the 
mean-valued eigenvalue problem. This results in a 
great computational saving. 

When studying the effect of uncertainties of ran- 
dom variables on the fundamental natural frequen- 
cies, it is convenient to study their squared value, 
i.e., eigenvalues, Ang 16 gives a good description on 
statistical analysis. 

The mean value of the eigenvalue for random vari- 
ables is obtained by taking the expected value of Eq. 
(43) 

Mx = E[ A] = A° + ^E **/ Ele^) (50) 
»=l 7=1 

The variance of the eigenvalues is obtained as 

KarfA] = £[A 2 ] — (51) 

For symmetrically distributed independent ran- 
dom variables, 

Var[X] = (52) 

\ E E A " 

»=1 k — 1 


( 48 ) 

( 49 ) 


where 


E[s 2 A 

E[e 2 A\ 

E[el\E[ S l\ 


= E 


<7=1 

N na „ 


{Xq ~ fix) 2 
N,am ~ 1 


( Xg — jJ. x )~{Xg — 


<7=1 

m ^ *am 

E E 

<7= 1 r =l 


(Xq - Vx)(Xr - M k 




Nsam 1 


and N sam is the number of samples, equal to ten 
thousand in this work. 

The standard deviation is calculated as 


< t a = ^V^X] ( 53 ) 


Calculating derivatives 

The above formulation only requires the calcula- 
tion of the derivatives of the stiffness matrix. These 
derivatives are obtained by taking the derivatives 
of the equivalent bending-stiffness matrix, eq. (18). 
However, the derivative for the first two rows and 
columns (reduced bending-stiffness matrix) are more 
involved. Various numerical schemes exist to eval- 
uate these derivatives. When using some of these 
numerical schemes, ill-conditioning could be a prob- 
lem. This problem can be avoided by the following 
formulation which allows the derivatives to be ob- 
tained exactly by numerical multiplication. The 
technique consists in taking derivatives of eq. (17) 

(D«U - [0/./U - (54) 

[Di,u\ Xt [Dtijr] 1 [D/ij\ - 
[Di.ii] [Dii.ii] x [ 1 D u,i] ~ 

[Di.ii] [Diui }~ 1 [ E>u.i } a , 

UMx . x , = 1 D './ Lx , " < 55 ) 

[D[,n\ iXtXi {D lul r'[D,u\- 
[Di.u] iIm [£>//,//]“' [Du. /I - 

[Di.ii] x> [Duji\ 1 [Du.i] X) - 
[Di.ii] X) [Drui] [Dnj\ - 
[Di.u\ [Du.ii\Ax, Xl [ Du.i } - 
[Di.ii] [D , /.//]"! [Dtu\ Xi - 
[Di.u] Xi [Dii.ii] 1 [D/u] Xi - 
[Di.ii] [Dii.ii]~ X) [Du.i] It - 
[Di.ii] [Dii.ii] l [Du.i] XtX 
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Dimens*on!«ss Fundamental Eigenvalue for a 
9<V— 45/-45/90 Laminated Beam 


Dimensionless Fundamental Eigenvalue tor a 
9<y-45/0/25 Laminated Beam 




Figure 1: Dimensionless Eigenvalue distribution for \1CS 
and SFEA for a thin symmetrically 4-ply laminated 
beams 


Figure 2: Dimensionless Eigenvalue distribution for MCS 
and SFEA for a thin unsymmetrically 4-ply laminated 
beams 


The derivatives of [£//,//] 1 are calculated using the 
following matrix definition 

[Du,u] 1 [£>//,//] = (/] (56) 

and these derivatives being 

{Du,u]~ x \ = (57) 

~[Du,ii\ 1 [D/i'ii] Xt [Duji] 

I d ii,ii\~ x%Xj = (58) 

-[Du.n\~ x [ [Diiji] Xj [D//.//] -1 - 

[Du.ii]-^ [Duji] [Dii'ii] 1 - 

[Du.u] 1 [Dti'U\ ZXi [Dii'U\ 1 


Dimensionless Fundamental Eigenvalue tor a 
9Q/-4 Srawa/3<y-45tf0 Laminated Beam 



RESULTS 

The numerical results are obtained for a can- 
tilevered thin laminated beam. The plies are as- 
sumed to be made of Graphite- Epoxy. The beam’s 
material and geometrical properties used in the anal- 
ysis are§: 

Afateria/ density: 


P 


1.449 x 1(T 4 


slugs 

in 3 


Major in-plane Poisson’s ratio: 


Figure 3: Dimensionless Eigenvalue distribution for MCS 
and SFEA for a thin symmetrically 8-ply laminated 
beams 


Young’s Modulus: 

Eyy = 1.36 x 10 6 psi 

E xx = 13.75 Eyy 

Shear Modulus: 


u xy = 0.30 

5 All properties are assumed to be uniform throughout the 
beam 


G X y 

= 0.55 Eyy 

Gyz 

= 0.25 Eyy 

G xz 

= 0.25 Eyy 
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Table 1: Statistics analysis 



MA 

X 

b 

10“ 2 

Ply-orientations 

MSC 

SFEA 

MSC 

SFEA 

[90/-45], 

1.63468 

1.63925 

0.41519 

0.41928 

[90/-45/0/25] 

4.10664 

4.12107 

6.44410 

6.67425 

(90/-45/30/0), 

3.23110 

3.23559 

2.79163 

2.80235 

[90/-45/30/0/ 

/25/45/-90/-30] 

4.53781 

4.54736 

9.51499 

9.69171 


Dimensionless Fundamental Eigenvalue for a 
90/— 45/3CVQ/2S'45/-9G/-30 Laminated Beam 



Sensitivity o* the Fundamental Eigenvalue for a 
9CV-45/-4V90 Laminated Beam 



Sample number 


Figure 4: Dimensionless Eigenvalue distribution for MCS 
and SFEA for a thin unsymmetrically 8-ply laminated 
beams 

Beam's dimensions 

Width : b Q = 1.0 in 

Thickness : h Q = 0.5 b Q 

Length . — 30.0 /ig 

For both the stochastic finite element analysis 
and Monte Carlo Simulation ten beam elements 
were used and four different laminated beams are 
considered: (i) symmetrically laminated cantilever 
beam with four plies, (ii) and unsymmetrically lam- 


Figure 5: Sensitivity of the dimensionless eigenvalue for 
a thin symmetrically 4-ply laminated beams 

inated cantilever beam with four plies, (iii) symmet- 
rically laminated cantilever beam with eight plies, 
(iv) and unsymmetrically laminated cantilever beam 
with eight plies. 

The natural frequencies are non-dimensionalized 
with respect to the deterministic quantities as fol- 
lows 

- ujnLt * am nr (59) 

ho V Ey y 

In general, the ply-angle uncertainties are between 
±2°. However, results show the effect in considering 
twice the uncertainty in ply orientations, i.e., from 
-5° to 5°. The statistical analysis is shown in Table 
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SeostOvity of th« Fundamental Eigenvalue for a 
9Q/-45/Q/25 Laminated Beam 



Sampie number 

Figure 6: Sensitivity of the dimensionless eigenvalue for 
a thin unsymmetrically 4-ply laminated beams 


Sensitivity of frie Fundamental Eigenvalue for a 
90/-45/30/0rtV3(V-i5«0 Laminated Beam 



-05 


o 1000 2000 3000 4000 5000 6000 '000 *000 9000 0000 

Sample numoer 

Figure 7: Sensitivity of the dimensionless eigenvalue for 
a thin symmetrically 8-ply laminated beams 


1 . 

The results in Figures (l), (2), (3), and (4) show 
that the present model has a good correlation when 
compared to MCS. These results are obtained for 
only 1000 samples opposed to ten thousand us- 
ing MCS. As the number of plies is increased the 
model correlated even better with MCS. The mean- 
valued fundamental eigenvalues are accurately ob- 
tained with one hundred samples using the present 
model opposed to ten thousand Monte Carlo Simu- 
lations. 

Figures (5), (6), (7), and (8) show that the eigen- 


Seneitrvity of me Fundamental Eigenvalue for a 
9<y-45/3CMy25/45/-9Q/-30 Laminated Beam 



Semple number 

Figure 8: Sensitivity of the dimensionless eigenvalue for 
a thin unsymmetrically 8-ply laminated beams 


values are sensitive to the first derivatives and the 
second derivatives are not influential. These fig- 
ures show that the sensitivity analysis is significant, 
therefore the ply angle uncertainties can play an 
important roll in affecting free vibrations of symmet- 
rically and unsymmetrically laminated thin beams. 

CONCLUSIONS 

Monte Carlo Simulation has been applied to thin, 
symmetric and unsymmetrically laminated beams 
with randomness in ply orientation to study the free 
vibrations. At least ten thousand realizations of the 
Monte Carlo sampling have been performed to im- 
prove the accuracy of the analysis. 

A second stochastic finite element approach has 
been developed using perturbation methods. Using 
Taylor Series expansion the eigenvalues has been ex- 
pressed as mean-valued and probabilistic quantities. 
The accuracy of the results have been compared to 
those obtained by Monte Carlo Simulation. 

An elegant way to obtain sensitivity derivatives 
is detailed. The present method is advantageous 
over other techniques because the eigenvalue prob- 
lem needs to be solved only once. With only one 
hundred samples our model agrees with ten thou- 
sand MCS. 

Based upon the results, this method results in a 
great computational saving when one is interested 
in predicting the statistics of the fundamental natu- 
ral frequency of unsymmetric laminated beam in the 
presence of ply-angle uncertainties. 
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